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Abstract 

We consider the characterization of the nonequihbrium stationary state of a randomly- 
driven granular gas in terms of an entropy-production based variational formulation. 
Enforcing spatial homogeneity, we first consider the temporal stability of the sta- 
tionary state reached after a transient. In connection, two heuristic albeit physically 
motivated candidates for the non-equilibrium entropy production are put forward. 
It turns out that none of them displays an extremum for the stationary velocity 
distribution selected by the dynamics. Finally, the relevance of the relative Kullbach 
entropy is discussed. 
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1 Introduction 



Apart from being the subject of intense experimental activity, granular gases 
are also a particularly fertile testing ground for new theoretical ideas and 
problems, especially within the field of nonequilibrium statistical physics. One 
such a problem is the role of entropy production as a Lyapunov functional 
for nonequilibrium steady-states. This problem has its roots in the fifties, in 
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the works of the Brussels group around Prigogine [1] on the minimum entropy 
production theorem. The limitations of this theorem, that relies essentially on 
the linear response formalism (i.e., has a domain of validity that is restricted 
to close-to-equilibrium situations), were rather clear already at that time, and 
a first extension to far-from-equilibrium situations was proposed under the 
form of the phenomenological "general evolution criterion" of Glansdorff and 
Prigogine (see [2] and references therein). 

Since then, there was steady effort, and a huge body of literature, in the search 
for a variational principle for steady-states that are arbitrarily far from equi- 
librium. Several candidates for a "nonequilibrium entropy production" with 
extremal properties at stationarity were proposed in different contexts, and 
at various levels of coarse-graining of the description - from the microscopic 
to the phenomenological ones. Some success was encountered for Markovian 
systems described by a master equation for the probability distribution func- 
tion of the microstates - starting with the pioneering work of Jiu-li et al [3], 
and intensively studied afterwards (see, e.g., [4,5] to cite only a few). Also, a 
connection between the phase space contraction rate in dissipative, externally 
driven systems and an entropy production rate was established in some cases, 
see e.g. [6] for a critical discussion. An extension of Jaynes' maximum entropy 
inference principle (MaxEnt) to nonequilibrium situations was proposed [7], 
and illustrated recently on several examples [8]. The resulting picture is, how- 
ever, rather confusing and sometimes even contradictory (e.g., some of the 
above-mentioned papers speak of a "maximum" entropy production rate at 
stationarity, while others refer to a "minimum"). 

One of the main difficulties of nonequilibrium statistical mechanics is the 
scarcity of solvable models, on the basis of which one could, eventually, get 
some clarification on these controversial points. The purpose of the present 
work is to consider such a solvable model, namely a granular gas modeled as 
an assembly of inelastic hard-spheres with constant restitution coefficient, in 
which energy is injected by means of random forces acting independently upon 
the particles. The balance between dissipation and the random kicks allows 
the system to reach a nonequilibrium steady-state (NESS). In a Boltzmann 
equation description, one can compute (in some perturbative expansion) the 
single-particle probability distribution function (pdf). This model is widely- 
used and very successful in explaining many features of granular systems (see, 
e.g., [9,10,11]). One of the question is thus whether this model is also ap- 
propriate in describing thermodynamical properties of granular systems - in 
particular, the entropy production rate and its eventual relationship with the 
relaxation to NESS. We propose two heuristic - albeit physically motivated - 
candidates for the nonequilibrium entropy production rate, as functionals of 
the pdf, and we discuss their extremal properties in NESS. Such a granular 
gas has a strong "built-in" irreversible element at the very level of the grain 
dynamics, which is represented by the inelasticity of the collisions. However, 
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one may ask whether in the hmit of a very weak inelasticity (i.e., for steady- 
states that are arbitrarily "close to equilibrium") one could recover a kind 
of "minimum entropy production theorem" in a stochastic formulation - an 
equivalent of that described in [3]. We will also address this point here. 

In the next section we are introducing the model, and in Sec. 3 we study the 
nonequilibrium steady-state and its linear stability. Section 4 is devoted to the 
discussion of the nonequilibrium entropy production issue, and the behavior 
of the relative Kullback entropy. We conclude in Sec. 5 with a brief discussion 
of the limitations of this model as far as describing the thermodynamics. 



2 The model 

We consider a granular gas of inelastic hard spheres in dimension d ^ 2, 
uniformly heated by a stochastic thermostat, as described in detail in [9,10]. 
The particles undergo binary inelastic collisions, modeled through a constant 
restitution coefficient a G [0, 1] that is meant to characterize the degree of 
inelasticity; the limit a = 1 corresponds to perfectly elastic collisions, while 
a — corresponds to the perfect inelastic ones. Each particle i (of mass m) 
is subjected to an external Gaussian white noise force ^j(i); these forces are 
uncorrelated for different particles, and homogeneous in space, 

{^iAtKjAt')) = m^^^oS^ApS{t -t'), a,P=l,...,d. (1) 

We describe the system at the level of the kinetic theory, and for simplicity, 
without affecting the overall conclusions, we shall concentrate on the spatially 
homogeneous case. For the single particle distribution function /(r, vi, t) ~ 
/(vi, t), the Boltzmann equation reads then: 

9*/(vi,t)=x/[/,/] + f^/(vi,t). (2) 

The extra term (^Q/2)(9^/9vf)/(vi, t) accounts for the change in the distribu- 
tion function caused by the random "kicks" the external thermostat is applying 
on the grains. It corresponds to an injection of energy at constant rate 
per unit mass, x is the pair correlation function at contact and 

nf, /] = / / d^ e{a ■ vi2)(^ • V12) 

X (lr^-l)/(vi,t)/(v2,t) (3) 

represents the inelastic two-particle collision operator. Here a is the diameter 
of the hard spheres; a is a unit vector joining the centers of the particles at 
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contact; V12 = vi — V2: 0{...) is the Heaviside step-function; and 6 is an 
operator that restitutes the pre-collisional velocities, i.e., 



1 + a 

6~Vi = vf=Vi — (vi2-5-)o-, (4a) 

Za 



b V2 = V2* = V2 + ^ J" (vi2 • CT)a- . (4b) 

za 

Note that the post-collisional velocities are 



6vi = V* = Vi - ^-^(vi2 •a-)a-, (5a) 
= V2 = V2 + ^ "1 (vi2 • a-)a- . (5b) 



3 Scaling solution and stationary state 

3. 1 Scaling solution of Boltzmann's equation 

It turns convenient to introduce the pdf / of rescaled velocities c — v/vt- 
where n is the number particle density and 



-W^^ (T, 
is the thermal velocity associated to the kinetic temperature of the particles, 

\ksnt)^^- / dv%V(v,t) (8) 
{Ub is Boltzmann's constant). 

For inelastic collisions, /(c, t) is different from a Gaussian 
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and it is customary to characterize its deviation from a Gaussian through a 
series development in terms of Sonine polynomials 5'„(c^), which, in practice. 
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(10) 



(11) 



is truncated to the first non-zero term [11], 

f{c,t) = <P{c) [l + a2{t)S2{c') 

where 

2 2 8 

The possible explicit temporal dependence of /(c, t) appears through the time- 
dependent coefficient a2{t) of the Sonine polynomial 5'2(c^). 

For consistency of the description, it is foimd that the kinetic temperature 
T{t) and the coefficient 02 (t) obey a set of two coupled nonlinear first-order 
differential equations: 



dT{t) kBnxa''-\l-a'')Sd 



dt 



7rm 



d 



X 



1 + 



1024 



al{t) 



, (12) 



da2(t) 2mil 



dt 



kBT{t) 



X 



a2{t) + 



^kBT{t) Anxa''-\l-a')Sd 



TTTJl 



d{d + 2) 



—aJt) H —an(t) 

16 ^ ^ 1024 2^ \ 



nm d{d + 2) 



1 + 



' d{d + 2) 

1 + g Q2(^) 



+ Dia2{t) + D2al{t) 



(13) 



Here Sd = 2%'^^'^ /T{d/2) is the surface of the unit-radius sphere in dimension 
d, r being Euler's Gamma function. Equation (12) follows from the defini- 
tion (8) of the kinetic temperature, while Eq. (13) is obtained from the limit 
of vanishing velocities of the Boltzmann equation (2), see [12]. The coefficients 
Di and D2 are given, respectively, by [12]: 



l-2d-d'^ 



2{l+ay(d'-2d-5) 



8(l+a2)3 

+4((i-l)(a-l)2(l + a2) + 8(a^ + 6a2 + i)l , (14) 



Do 



d{d+2) 



+ 



12a^{l + a^){d-l){d-2) 



64 ' 32(l + a2)5 

-Aa'^{l+a^){2A+Ad-d'^) +Aa{l+a^){d+6){d-l) 

-{l+a^){26+28d+9d'^ 



(15) 
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3.2 Stationary state 



In the asymptotic limit, the granular system will reach a stationary state, 
that results from the balance between the energy injection by the external 
thermostat, and the energy dissipation through inelastic collisions between the 
particles. The stationary temperature Tq is thus related both to the restitution 
coefficient a and to the amplitude of the Gaussian thermostat. Or, to state 
it differently, in order to ensure a given value of Tq (for a fixed value of a) , as 
resulting from the stationary form of Eq. (12), one has to tune the amplitude 
of the stochastic thermostat to 



nxo- 



d-li 



1 - a^)Sd ( ksTo 



3/2 



m 



X 



' + 16^^° + 1024'^^'^ 



(16) 



Here 020 is the stationary value of the coefficient of the first correction to the 
Gaussian. Its expression can be obtained from the stationary form of Eq (13) 
and it is the solution of the third order nonlinear equation (see e.g. [11] for 
a discussion concerning the relevance of the corresponding three roots in the 
case of a force-free system): 



:i-«') (1 + ^020 + 
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1024 



^20 



1+020 



(d+2)(d+4) 



= V2 



-a 



+ Dl 020 + i^2 020 • (17) 



The coefficient 020 can be obtained in a closed analytical form through a 
Taylor expansion of the above equation. It was however shown in previous 
works [12,10] that there are some ambiguities from this linearization procedure 
that may affect 020. We therefore chose the linearizing scheme that yields the 
closest result to the Monte Carlo simulations of Ref. [12]: 

= -16(1 - q;2)(1 + a^){l-V2 + a^) 

X {l6V2+13+4d(3V^+l)+2d^(V^-l) 
+o;^(— 75+44d— 2d^) — CK^ 

+Q;^(-5 + 4d+2d^)}"\ (18) 

Considering instead the expression derived by van \oije and Ernst in [9] would 
not alter the following discussion. Note that 0,20 becomes zero in the elastic 
limit a = 1, when the stationary probability distribution recovers trivially the 
Gaussian, equilibrium shape. 
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The corresponding stationary probability distribution function is therefore 



n 



n 



/o(v) = ^/o(c) = ^0(c) 1 + a2o^2(c ) 



'To 



'To 



(19) 



where vtq = \/2kBTo/m is the stationary value of the thermal velocity. 



3. 3 Linear stability analysis of the stationary state 



The stability of the aforementioned steady state has not been investigated in 
the literature, even if the hydrodynamic-like equations have been derived re- 
cently for the (dilute) system considered here [13]. A complete linear stability 
analysis (and its eventual comparison with the existing results for the homo- 
geneous cooling state [14]) is a tedious task, and a separate research subject 
that we shall not address here further. Instead, we shall consider a simplified 
version of it, in which the homogeneity of the state is not affected by the 
perturbations. This will by no mean influence our general conclusions. 

Let us then consider small deviations of the temperature and of the coefficient 
02 from their stationary values, 



T = To(l + 59) , 02 = 020 + 5a2 , 
with \56\ <^ 1 , \Sa2\ <^ \0'2o\- 



(20) 



The linearized evolution equations of these perturbations result from Eqs. (12) 
and (13), 



3 3/16 + (9/512)o2o 

2 l + (3/16)o2o + (9/1024)oio 



(502 



(21) 



020 59 



d + 4 



1 + 020 



d{d + 2) 



V2 



l — oi' 



(L>i+2L»2a2o) 



3 9 

1 ^20 

16 512 
.3 2 

H 020 H Onn 

16 1024 



502 



f22) 



In Fig. 1 we have represented the two eigenvalues of the corresponding stability 
matrix as a function of the restitution coefficient a, for both d = 2 and 0? = 3 



cases. 
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Fig. 1. The eigenvalues of the hnear stabihty matrix for the stationary state as 
a function of a for d = 2 and d = 3. The eigenvalues are measured in units 

One notices that the two eigenvalues are strictly negative for a < 1, which in- 
dicates the stability of the stationary state with respect to small perturbations 
in the temperature and in the shape of the probability distribution function 
(in the scaling form) . We emphasize again that spatial homogeneity has been 
enforced here. As expected, in the elastic limit a = 1 one of the eigenval- 
ues becomes zero (while the other one remains negative) - which corresponds 
to the temperature becoming a marginal mode, and to a relaxation of the 
distribution function to its equilibrium Gaussian shape. 



4 Entropy production 

We now turn to the issue of entropy. For our homogeneous system, we consider 
the Shannon information entropy 

5(t)^-A;,/^^dv,/(v„^)ln(:^) (23) 

(with Euler's number e and Planck's constant /i, and eh'^ the volume of the 
semiclassical elementary phase-space cell)." It is known that in the elastic 
limit a — 1 (and in the absence of an external drive) this reduces to the 
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appropriate expression of the usual thermodynamic entropy and leads to the 
classical "//-theorem". We now consider the time evolution of S{t) as governed 
by the Boltzmann equation (2), which reads 



dS_ 
'dt 




e¥ J 



] 



(24) 



As mentioned in the Introduction, we wish to introduce a heuristic - albeit 
physically motivated - entropy production functional that, hopefully, displays 
extremal properties in NESS. We shall propose two approaches. But before 
proceeding further, we would like to remind the reader the status of dS/dt 
within the framework of phenomenological thermodynamics as discussed in 
standard textbooks [15,16,17], as well as some of its extensions to stochastic 
systems [3,4,5]. Entropy variations are usually split into two parts: 



where a-^-r ^ is the entropy production arising due to the dissipative pro- 
cesses that take place inside the system (that is positively-defined according 
to the second principle of thermodynamics), while the entropy fiux (Jfjux = 
~ Iv dW ■ Js accounts for the external forces driving the system into a 
nonequilibrium state (the related contribution is often reduced to boundary 
terms). The "art" of phenomenological thermodynamics precisely bears on Js 
and on how to decompose it in terms of the energy, particle, momentum, 
chemical, etc., currents. This is done, usually, on the basis of the local equi- 
librium hypothesis. In a similar way, ai„ often appears as a bilinear form in 
the fluxes rimning through the system and the conjugate affinities. In the 
near-to-equilibrium regime, the fluxes are usually proportional to the conju- 
gated affinities, with the Onsager coefficients as proportionality factors, and 
one recovers Prigogine's minimum theorem for ai„ under the hypothesis of 
time-reversibility of the underlying microscopic dynamics. 

However, in view of the local character of the energy injection mechanism, as 
well as of the spatial homogeneity of the system, the situation is completely 
different in the case we are considering. Indeed, unlike the above-mentioned 
"conventional" NESS, there are neither macroscopic, however weak, currents 
running across the system, nor the related phenomenological Onsager response 
coefficients. Therefore, the separation into "source" and "flow" for the entropy 
variation is much more tricky. 



dS _ 



(25) 
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First approach. - A first proposed choice of the "entropy production" is 



C^irr = j dVidV2da-^(o- • Vi2)o- • Vi2 



X (/r/2* - /1/2) In (^^] + I / dv^^^ , (26) 



the form of the first r.h.s. term being simply chosen by analogy with the elastic- 
limit case. The second term has been chosen by analogy with standard diffusion 
processes. In those processes this term vanishes at equilibrium because the 
gradients disappear. Note however that here the diffusion process happens 
in the velocity space, and thus the vanishing of this term at equilibrium is 
not due to the system becoming homogeneous in v-space, but because of the 
energy source strength being tuned to 0. The above (Ti„ appears to be the 
sum of two positive definite terms, and it is therefore also positive definite. 
Furthermore, cxi^^ can only be zero at equilibrium, namely when both the energy 
source (the random kicks) and the energy sink (the dissipative collisions) are 
tuned to zero. In that respect, it fulfills the properties expected from standard 
phenomenological thermodynamics. 

On the other hand, the form of the entropy flux CTflux is now constrained to be 



efflux = — J dvidv2d6-6'(o- ■ ^12)0- ■ V12 

"(/r/2**)(/i/2) 



X /1/2 In 



*\2- 



(27) 



where we have used the shorthand notations /i,2 = /(vi^2,^), respectively 
/r,2 = /('^i*2) ^) fo'^' distribution functions corresponding to the pre-coUisional 
velocities (4) . The above functional of / is negative for a large class of trial 
functions, and must definitely assume a negative value Uflux ~ —^—^Tq^"^ 
in the steady state {i ~ ^Ji-i is the mean free path). However, aside from 
conveying the shrinking of phase space volumes, we must dismiss (Tirr/flux as 
relevant candidates for extremum entropy functionals. Indeed, in the spirit of 
phenomenological thermodynamics, the splitting of dS/dt into ai„ and ctaux is 
motivated by the desire to isolate the driving processes (the source and sink 
referred to above) from the irreversible processes inside the system. However 
there is no simple and univoque manner to do so, and definitely this first 
choice is not accomplishing this physically-motivated requirement. It must be 
noted that the last term of eq.(26) could have also chosen as a part of cranx, 
which would then have featured both the source and the sink, at the price of 
abandoning its negative definiteness. 

Second approach. -We now propose an alternative and perhaps more pragmatic 
route, which consists in isolating as the only driving mechanism the random 
kicks provided by the thermostat. The inelastic collisions, viewed above as an 
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energy sink, are now incorporated into a term describing the system's intrinsic 
dissipative microscopic dynamics. Along those lines we henceforth write that 



dS . 

-TT — Cr syst + CText ■ (28) 

at 

The first contribution asyst corresponds to the entropy production inside the 
system, i.e., it comes from the changes of the particles velocities during the 
binary inelastic collisions. 



d-i 

syst 



a 



/ dvi / dv2 / dd-9{d- ■ V12) 



X {a- ■ Vi2)/i/2 In (^-H j , (29) 

where we have used the shorthand notation — f{'^i,2j't) for the distribu- 
tion functions corresponding to the post-collisional velocities (5). Of course, in 
the limit of elastic collisions a = l the expression of o-gyst reduces to the usual 
positive-definite expression of the hard-disk gas that enters the if-theorem. 
However, in general asyst does not have a definite sign. One can imagine the en- 
tropy production inside the system as resulting from two antagonist (although 
actually undissociated) mechanisms, namely a generic disordering effect of any 
particle collisions (e.g., that is also present for elastic hard spheres) in d ^ 2, 
and an ordering effect due to the inelastic character of the collisions (i.e., to 
the reduction of the translational agitation of the particles). Depending on the 
actual shape of the distribution function, one of these two mechanisms may 
prevail on the other, thus determining the sign of the instantaneous value of 

^syst ■ 

The second contribution (Text is determined by the effect of the thermostat on 
the distribution function of the particles of the system. It corresponds to an 
energy injection into the system, and to a disordering effect of the particles 
velocities (through "random kicking"), and therefore, as expected, is always a 
positively-defined quantity. 



Introducing the dimensionless quantities 

~ "^^systjext 

(Tsyst,ext = ^^d-l^^y ' ^^^) 

one obtains the expressions for the dimensionless time-dependent entropy pro- 
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duction sources: 



' syst 



m 

To 



1/2 



/ dci / dc2 / dd-6(a- ■ C12) 
X (o--Ci2)/(ci,i)/(c2,i) In 



/(Ci,t)/(c2,t) 

[f{cl,t)f{clt)\ 



respectively 



(32) 



To 
T{t) 



2d\/2n 



3 9 

IH a2(i)H al(t) 

16 ^ ^ 1024 ^ 

X / dc 

^1 



H'^ f{c,t) 



Ve/(c,t) . (33) 



In the stationary regime at temperature Tq one has, obviously, agyst = —^ext = 
— ctq. The quantity ao is positive and decaying monotonously with a, as illus- 
trated in Fig. 2. Note that ao is nonzero as long as the collisions are inelastic, 

6 



;b 3 
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0.6 
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a 



Fig. 2. The dimensionless stationary entropy production ctq as a function of a in 
d = 2 and d = 3. 

i.e., as long as the stationary probability distribution is non-Gaussian. Note 
also the negativity of (fsyst in the stationary state - the ordering effect due 
to the inelastic character of the collisions prevails on the generic disordering 
effect of the collisions. 
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Let us now address the question whether the entropy production (as a whole, 
or one of its parts (fgyst or o'ext) can play the role of some kind of "nonequi- 
librium potential" for the system, i.e., whether or not it can account for the 
linear stability of the stationary state of the system. The particular case of 
the quasi-elastic limit e = l — a<^lis especially interesting, given that 
the stationary state is close to equilibrium. One might then expect a priori 
that a "minimum entropy production theorem" (in the spirit of the "extended 
Prigogine theory" [3]) might be valid in this case. 

Consider thus small perturbations of the temperature and of the coefficient 02 
around their stationary values, as in Eq. (20). A Taylor development of the 
entropy production terms o-gyst and aext leads to nonzero linear contributions 
in the perturbations 59 and 5a2, 



'syst 



-ao)^-S9 (^)+<5a2/^/c,/^^dc2/da 

S2{cl) 



X e{a- ■ Ci2)(o- ■ C12) /o(ci)/o(c2) 



X 



+ 



S2{4) S2{cf) 



1 + 020.^2(01) l + a2oS2{cf) l + a2oS2{c*2 



1 + 020.^2(0?) 

,*2^ 



S2{4) 



1 + 020-52 (cf) 1 + 020-52(0!) 



In 



/o(Cl)/o(C2) 

L/o(ct)/o(c^) 

+ 0(66^64), 



(34) 



respectively 



cText - {ctq) = -S9{ao) + Sa2< 



3/16 + (9/512)020 
1 + (3/16)020 + (9/1024)aio''° 



+ 



:i-a^)Sd . 3 9 



H 020 + 

2c^V27r V 16 '° 



1024^2° 



X / dc 



2 (Vc/o(c)) ■ [Vc ( 



e-'S2(c^ 



0{5e^,6al). (35) 



The total entropy production (Tgyst + f^ext also contains linear terms in the 
perturbations 56 and 5a2. 

The same holds true even in the quasielastic limit £ = 1 — a -C 1, when 
one can evaluate explicitly to O(e^) the expression of the coefficients of the 
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perturbations. More precisely, 



(-5^0) = -59 (^) - 5a, 



V{d/2) 



X 2{d - 1) + 



= -5e{do) + 5a2 ^^o + 0{e'') , (36) 




O'syst — Co 



where the stationary values are 



e + 0(£') , 



(37) 



and 



020 = 



\/2(V2-l) 
d-1 



e + 0(£2) . 



(38) 



The meaning of this result is that the entropy production as defined above 
cannot be used for a variational description of the relaxation of the system 
towards the stationary state, not even in the quasi-elastic limit. 

One may argue that the choice of the definition of the entropy production 
inside the system might be inappropriate, since it refers only to the trans- 
lational degrees of freedom, and it does not take into account properly the 
internal degrees of freedom of the particles - that are, in fact, responsible 
for the inelastic character of the collisions. The description of the inelasticity 
through a constant restitution coefficient a might thus be incompatible with 
a thermodynamic description of the system in terms of entropy production. 
We note that it is known that such a model, although being a useful approxi- 
mation which captures important physical effects, is in fact incompatible with 
basic mechanical laws (see e.g. chapter 3 of Ref . [11]). 

Let us now discuss briefly another issue that draw attention recently, see 
Refs. [18,19], namely that of the KuUback relative entropy, defined as 



Sji{t) is a measure of the "distance" between the actual pdf /(v.t) and its 
stationary profile /o(v), and, of course, is equal to zero at the stationary 
state. Following Ref. [19], one can parametrize SR{t) through the two sets of 
parameters, {71 = a2(t),72 = T{t)} for the nonstationary state, respectively 




(39) 
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{7io = O20 - 720 = ^0} for the sationary pdf. Considering as above (sec. 3.3) 
small deviations of the temperature and of the coefficient 02 (that result in a 
small deviation 6f{v, t) of the pdf ) from their stationary values, one finds: 

SSn^-^ [ .^^J^ASfM]' = E Fij5^^S^2 ^0, (40) 

2 Jwi /o(v) ^i,j=l,2 

where 5ji — 5a2, 6^2 — Tq59, and Fij is the positively-defined Fisher informa- 
tion matrix [20], 




It looks therefore as if this relative entropy has the required property of ex- 
tremum at the steady-state (and monotonous exponential asymptotic relax- 
ation towards it). This property has already been demonstrated for other types 
of nonequilibrium stochastic systems (e.g., in Ref. [18], the one-dimensional 
Ornstein-Uhlenbeck and Rayleigh processes, noise-perturbed harmonic oscil- 
lator, dichotomous noise). The question arises about its relationship with the 
thermodynamic entropy production; in Ref. [19] it was shown that in the 
case of the usual Smolukhowski diffusion the Kullback entropy time- variation 
rate coincides with the Shannon entropy production rate. However, some fur- 
ther case-study (in particular, on systems described by kinetic Boltzmann- 
like equations) are necessary before generalizing this important conclusion to 
other nonequilibrium situations. In particular, although very appealing, the 
Kullback entropy does not reduce to the usual i?-functional in the limit of 
an elastic gas of particles relaxing to equilibrium. Besides that, computing 
Kullback entropy requires the knowledge of the steady-state pdf, while the 
expected approach would be to define a proper Lyapunov functional of the 
system from which to deduce the stationary state. 



5 Conclusions 

We illustrated on the well-known model of a randomly driven granular gas 
with constant restitution coefficient the difficulties that one encounters when 
trying to construct a variational principle for NESS based on an "entropy pro- 
duction". Two approaches were proposed for the interpretation of the entropy 
balance equation in terms of "sources" and "flows", but none of them lead to 
the formulation of such a principle. The main reason for this failure seems to 
be the intrinsic irreversible microscopic dynamics of the granular gas. Model- 
ing the internal degrees of freedom of the grains (that are responsible for the 
inelasticity of the collisions) through a constant restitution coefficient is thus 
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not appropriate for a thermodynamic description. This shows thus a major 
limitation of this model. A further step in the rather involved question of the 
refinement of the description would be the use of random restitution coeffi- 
cients (as done, e.g., in [21]). These are meant to describe the possible flow of 
energy (at the collision) both towards and from the internal degrees of free- 
dom to the translational degrees of freedom. Such a model, however, cannot 
be treated analytically, and no simple analytic conclusions can be therefore 
drawn on the fate of the corresponding if- functional. Numerical results are 
left for further studies. 

Moreover, the problem of the Kullback relative entropy, its monotonous relax- 
ation to the steady-state, and its relationship with the thermodynamic entropy 
production of a nonequilibrium system is a very promising direction for further 
studies. 

We thank Ph. Martin and J. Piasecki for suggestions and discussions. I.B., 
F.C., and M.D. acknowledge partial support from the Swiss National Science 
Foundation. 
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